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Abstract. Spectacular results of pulsar observations in 
X-ray and gamma-ray domains gave a new boost to the- 
oretical models of pulsar magnetospheric activity. A chal- 
lenging aspect of these efforts is that lightcurves and 
broadband energy spectra of the brightest HE sources ex- 
hibit unexpected richness of features, making each object 
almost a unique case to be interpreted with a custom- 
made model. This review offers our subjective account of 
how the models of high-energy radiation used within the 
framework of SCLF polar-cap scenarios tackle with these 
challenges. We describe major characteristics of all radia- 
tive processes relevant for modeling the observed features. 
Then we address successes as well as noticeable disadvan- 
tages of these models upon their confrontation with the 
available data. 



1. Introduction 

A large fraction of galactic X-ray and gamma-ray sources 
is associated with neutron stars, in particular with rota- 
tion powered pulsars (RPP). Spectacular results of ob- 
servational campaigns of RPP with ROSAT, ASCA and 
CGRO, and recently with RXTE and CHANDRA induced 
a new wave of interest in theoretical models of pulsar mag- 
netospheric activity. A challenging aspect of these efforts 
is that lightcurves and broadband energy spectra of the 
brightest HE sources exhibit unexpected richness of fea- 
tures making thus each object a unique case pending spe- 
cial treatment. 

The pair creation paradigm is a pivotal element in vir- 
tually all models of magnetospheric activity of radiopul- 
sars (which at the moment constitute a vast majority 
of known RPPs). Electron-positron pairs (e -pairs) are 
thought to be responsible for radio emission observed in 
radiopulsars which is interpreted as the coherent curvature 
radiation of plasma. Pairs can be produced in magne- 
tospheric environments either via photon absorption in a 
dense field of soft photons (photon-photon collision) or 
via photon absorption in a strong magnetic field. In ei- 
ther case a supply of hard gamma photons is required 
in order to fulfill stringent threshold conditions for pair 



creation. Not all of these photons would be subject to 
absorption; many will escape the magnetosphere without 
any attenuation. This leads us to expect that all radiopul- 
sars (including law-B millisecond objects) should be the 
sources of gamma radiation. It is up to theoretical mod- 
els to show whether the expected flux of the radiation 
is interestingly high with respect to the sensitivity of re- 
cent and future gamma-ray telescopes. One should keep in 
mind, however, that the pair creation paradigm remains 
just a paradigm for the time being: recently it was ques- 
tioned by lessner et al. 2001, who presented arguments 



that for a wide range of surface temperature, magnetic 
field strength and spin period, the electrons supplied by 
neutron star surface via thermionic and field emission will 
screen out the accelerating electric field £u , limiting it to 
a residual value. In consequence, no favourable conditions 
would exist for magnetic pair production. 

This review offers our subjective account of how the 
models of high-energy radiation do perform within the 
framework of polar-cap models with space-charge-limited- 
flow (SCLF). First, we introduce all ingredients of polar- 
cap models with SCLF and we describe major character- 
istics of all radiative processes relevant for physical condi- 
tions in pulsar magnetospheres. Then we address successes 
as well as noticeable disadvantages of these models upon 
their confrontation with the data. 

2. Basic parameters 

For the sake of completeness, we recall in this section basic 
quantities used throughout the review. A neutron star of 
radius R s and moment of inertia I, spinning with angular 
velocity fl = 2tt/P which decreases in time (for whatever 
the reason) at a rate fl — —2ttP~ 2 P < 0, loses its rota- 
tional energy at the rate 
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10 31 7 45 A 



15 



P erg s 



(1) 



where P is in seconds, P_i5 = P/10 and Z45 = 
J/10 45 g cm 2 , and L sc j is the so called spin-down lumi- 
nosity. If a magnetic dipole is anchored at the center of 
a neutron star, with its magnetic moment inclined at 
angle a to the spin axis f2, then the strength of the field 
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at the stellar equator is B s = [ib/RI (and twice as high 
at the poles). The dipole, rotating in a vacuum will lose 
energy at the rate 



3c 3 



sin 2 a flf O 4 . 



(2) 



The quantity B s sin a can be thus inferred from P and P 
for a neutron star with known values of I and R s . Another 
model, where the dipolar radiation is replaced with a mag- 
netospheric wind of particles ( Goldreich fc Julian 19691 ), 
gives a similar result as (0) for an orthogonal case: 



L 



wind 



-4 Bt Rl n 4 



(3) 



and therefore is independent of the angle a. Since no ob- 
servational support exists for P depending on sin a, the 
standard approach is to apply the latter model (L s d = 
-^wind) to derive the strength of the dipolar component of 
magnetic field 



10 15 / 45 R^PP 



(4) 



where B 12 = B S /10 12 G, and R 6 = i? s /10 6 cm. It is 
likely that neutron star magnetic fields contain high- 
order multipolcs which may dominate the dipolar com- 
ponent at the surface level. Their relative amplitudes as 
well as distribution remain, however, unknown. It will be 
assumed throughout the paper that the magnetic field 
is a static-like dipole, not distorted by rotational ef- 
fects or by presence of strong outflowing wind of parti- 
cles (the latter effect has been recently invoked to de- 
crease very hi gh values of B s infe rred from P and P 
for two SGRs ( |Harding et al. 1999b ). The field is there- 
fore approximated with axisymmetric field lines satisfying 
r sin -2 — i?dc in polar coordinates r and 9, with the 
dipole constant Rdc- A dipole constant for which a rigid 
rotation with the angular velocity reaches the speed- 
of- light limit (it occurs for i?dc = c/fi and this particu- 
lar value is denoted as R\ c ) determines the so called light 
cylinder of radius R\ c . All field lines which cross the light 
cylinder are then considered as open lines, and their foot- 
points on the stellar surface define two polar caps of radius 
R pc « R s • (R s / Ric) 1 / 2 , where the latter factor is the sine 
function of the magnetic colatitude 8 pc for the outer rim 
of the polar cap: 



sin 8„ 



(Rs/ i?k 



a/2 



(5) 



3. Properties of pulsars in X-rays and 
gamma-rays 

For more than 1500 pulsars known to date only about 
~ 40 positive detections in X-rays and no more than 10 de- 
tections in gamma-rays have been achieved ( jBecker 2002 . 
Kanbach 2002). CGRO provided seven high-confidence 



detections (of Seven Samurai) with three other cases 
classified as 'likely' detections. The gamma-ray sources 



were identified by virtue of flux pulsations with previ- 
ously known P and P. Crab and Vela are the only pul- 
sars seen by all three instruments of CGRO. No trace of 
pulsed signal in VHE range (300 GeV - 30 TeV) has been 



found so far for the gamma-ray pulsars (Sako et al. 2000, 
|Wcekcs 2000 ). However, strong steady VHE emission is 



associated with 3 out of 10 gamma-ray pulsars. Two pleri- 
onic sources of the steady VHE radiation - The Crab Neb- 
ula and the plerion around B 1706-44 - m ay serve as sta n- 
dard candles, with 'grade A' according to Weekes 200C . A 
third plerion - around the Vela pulsar - was given 'grade 
B' in the same ranking. All 10 gamma-ray pulsars are 
strong X-ray emitters. 

The positions of these HE pulsars are shown in the 
P — P diagram of Fig.[l] along with positions of about 
700 radio pulsars for which P values were available. A re- 
markable fact is that the location of X-ray sources does 
not correlate with the inferred strength of magnetic field 
B s ; at least not in a naively anticipated way that high- 
B objects would emit HE radiation, whereas low-B ob- 
jects would not. In particular, 10 millisecond pulsars - 
about thirty percent of all millisecond pulsars (the ob- 
jects with P < 0.01s and P < 10~ 17 , i.e. with low B val- 
ues: B s 10 9 G) known to date - have been detected as 
X-ray sources. So far, millisecond pulsars have eluded the 
detection in gamma rays (possibly with one exception, see 
below) and just upper limits were available for a handful of 
them from EGRET observations flNel et al. 19961) . In the 



1/2 
sd 



case of J0437-4715 the upper limit is interestingly tight 
- in disagreement with the empirical relation L 7 cx L 
(see Fig.|. 

Empirical evidence that pulsars must somehow tap 
their high-energy activity from rotation comes from a 
strong correlation between the success of detection in 
X and/or 7-rays and the position in the lists of targets 
ranked by spin-down flux values L s d/-D 2 , which is pre- 
sented in a graphical form in Fig. |^. 

Spectral analysis for pulsars detected with ROSAT 
PSPC (0.1 keV to 2.4 keV) shows that in most cases a 
power -law spectral model provi des acceptable fits to the 
data ( Becker fc Truemper 1997 ). Moreover, an intriguing 
empirical relation between inferred X-ray luminosity and 
spin down luminosity was found, Lx ~ 0.001 L s d, confirm- 
ing rotational origin of most of the X-ray activity. An in- 
teresting point is that the relation was obtained for all the 
sources regardless their temporal characteristics (about 
50 % of all pulsars detected with ROSAT are unpulsed 
sources). A complementary empirical relation was found 
for pulsed emission from 19 pulsars observed with ASCA 
(0.6keV to lOkeV). Assuming opening angle of X-rays to 
be one steradian, the inferred pulsed X-ray luminosity cor- 
relates with spin-down luminosity as 

3/2 



L x = 10 



.34 



10 38 erg s" 



erg s 



(6) 



according to Baito et al. 1998 
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Fig. 1. P — P diagram for Rotation Powered Pulsars. 
The pulsars detected exclusively in radio are indicated 
with dots; they are taken mostly from the data base 



of Taylor et al. 1995. Thirty five pulsars emitting X-rays 
are indicated with bullets. These include two objects re- 
cently di scovered with RXTE : J0537-6910 in SNR N157B 
in LMC ( Marshall ct al. 1998| ) is the fastest young pulsar 
known, spinning twice as fast as the Crab pulsar but with 
similar value of spin down lu minosity; J1846-0258 in SNR 
Kes-75 flGotthclf ct al. 2000| ), with P = 0.32s, the highest 
P among all RPP. Seven bullets in circles indicate seven 
gamma-ray pulsars of high confidence. Dashed lines corre- 
spond to constant values of spin down luminosity L s( j. The 
upper line (L s d «4x 10 38 erg s~ 
sar and J0537-6910, the lower one (L s d ~ 
includes J2144-3933 - the slowest (P — 8.5s) radio pul- 
sar detected so far ( Young et al. 1999 ). Dotted lines cor- 
respond to constant values of the dipolar component of the 
surface magnetic field as inferred from P and P through 
(|) with R 6 = 1. 

A similar power-law empirical relation holds for 
gamma-r ays (cf. Fig.p|), but with a different power-law in- 
dex (e.g. Thompson ct al. 1999): 



4 x 10 38 erg s includes the Crab pul- 
3 x 10 28 erg s" 1 ) 



1(T 



10 38 erg s" 



1/2 



erg s 



(7) 



Important conclusion from (g) and (0) is that neither 
Lx nor L 1 becomes a sizable fraction of L s & . The most effi- 



Fig. 2. The ranking of RPPs by their spin-down flux 
Lsd/D 2 . The units of abscissa and ordinate are chosen 
to keep the resemblance with the P — P-diagram: the 
lines of constant spin-down flux (two such lines - con- 
tinuous and dashed - are shown as examples) have the 
same slope with respect to the abscissa (+3) as the lines 
of constant L sc j in Fig. [I]; the lines of constant B s (dot- 
ted) have the same slope (—1) with respect to the ab- 
scissa in both figures; the range of values for coordi- 
nates is also identical. Crab and Vela sticking out well 
above the continuous line posses the highest spin-down 
fluxes among the RPPs. The continuous line is trained 
on B1951+32 with L sd /D 2 w 7 x 1(T 8 erg s" 1 cm" 2 ; the 
dashed line is for the millisecond X-ray pulsar J0751+1807 
with L s d/D 2 «2x 10~ 10 erg s _1 cm -2 . Note a low over- 
all position of the gamma-ray pulsar B 1055-52. (From 



Kluzniak & Rudak 1996 , with recent updates) 



cient conversion of spin-down luminosity into high-energy 
radiation is taking place for B1055-52 - the oldest pulsar 
among Seven Samurai - with L 7 ~ 0.1 L s d- 

High energy lightcurves differ significantly from those 
in radio. Their most striking feature are relatively long 
duty cycles as well as phase shifts in comparison to the ra- 
dio pulses. Only for the Crab pulsar the peaks in gamma- 
rays as well as in radio wavelengths occur at the same 
rotational phases. The light-curve shapes fall into two 
categories. The Crab pulsar, Vela and Geminga show 
two sharp peaks separated in phase by 0.4 — 0.5 and 
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Fig. 3. Gamma-ray luminosity versus spin-down luminos- 
ity for seven pulsars (rilled dots) detected with the CGRO 
instruments. Opening angle of one steradian was as- 
sumed for the gamma-ray emission. Open triangles are the 



EGRET upper limits after Ncl ct al. 1996 for 350 objects 



including seven millisecond pulsars. T he filled triangle in - 
dicates the upper limit for J0437-4715 ( |Fierro et al. 1995|) . 
Note, that most of the upper limits are well above the 
maximum possible value for L 7 set by L 7 = L s d (dotted 
line) . The dashed line marks the empirical relation derived 

1/2 

for the CGRO pulsars: L 1 oc L s ' d . Solid lines show evo- 
lutionary tracks for a classical pulsar with B s = 10 12 G 
(upper line) and a millisecond pulsar with B s — 10 9 G 
(lower line) according to the phenomenological model of 
Rudak k Dyks 1998[ 



connected by an interpeak bridge of considerable level. 
B1706-44 shows two peaks separated by 0.2 in phase, with 
some hints of a third component in between. Other pul- 
sars exhibit broad single pulses. Unknown opening an- 
gles for gamma-ray emission introduce a factor of uncer- 
tainty when inferring the gamma-ray luminosities. Broad 
peaks in gamma-ray pulses do not necessarily mean large 
opening angles for gamma-ray emission. Polar cap mod- 
els, which rely on purely dipolar magnetic fields postulate 
nearly aligned rotators, where inclination of magnetic axis 
to spin axis is comparable to the angular extent of the po- 
lar cap ( Daugherty fc Harding 1994 ). 



With a dicovery of sharply peaked pulsed X-ray 
emission in the fastest millisecond pulsar B1937+21 an 
apparently separate group of millisecond X-ray pulsars 
emerges, with its members - B1821-24 (Saito ct al. 1997), 
J0218+4232 (jKuipcr et al. 200q ), and B1937+21 



(Takahashi et al. 2000) being scaled-down versions 



of the Crab pulsar as far as sharp pulse profiles and hard 
power-law X-ray spectra are concerned. An astonishing 
common feature within the group is the same strength 



of the magnetic field estimated at the light cylinder and 
the fact that it matches the strength of the Crab pulsar 
magnetic field at the light cylinder. 

It is, however, lightcurves and broadband energy spec- 
tra extending from optical to X-rays and to gamma-rays, 



available for Seven Samurai (e.g. Thompson et al. 1999, 
Thompson 2001) that are particularly impressive and 
challenging for any models of pulsar activity. These ob- 



servations are reviewed in details by Kanbach 2002 (this 
volume). 

4. Unipolar induction 

Consider a neutron star of radius R s and surface mag- 
netic field B s as a perfect conductor. For the rotating star 
with its dipolar magnetic field immersed in a vacuum an 
external quadrupole electric field £ develops, with non- 
zero component along magnetic field lines at the surface. 
The corresponding elec trostatic pot ential $ in polar coor- 
dinates r and 9 reads ( Michel 1991 ) 



$(r, 



~ 3cos 2 9- 1) 



(8) 



where Q = ir B s i?g/(3 c P) is the quadrupole moment. 
Maximal electromotive force will then be induced between 
the pole of the star and the outer rim of the polar cap due 
to: 



where 



equator 



l 



Rpc 



equator 



= —BQRi 
2c 



(9) 



(10) 



denotes a (huge but unavailabale) potential drop between 
the pole and the equator. In terms of voltage, polar caps 
offer then AV pc «7x 10 12 B 12 P~ 2 Volt 



4-1. Potential drops in SCLF gaps 

Potential drop (||) sets a firm upper limit for all (model- 
dependent) realistic potential drops A$n, inc luding vac- 
uum gap model of [Ruderman fc Sutherland , as well as 
outer gap models of e.g. |Cheng et al. 1986 In fact, 



A$n < A$ 



pc 



(11) 



for most cases. We will concentrate hereafter on a class 
of models where the supply of charged particles from a 
neutron star surface along open field lines is not limited 
by cither binding or cohesive energy of the particles and 
therefore can reach the Goldreich-Julian rate close to the 
surface ('Space Charge Limited Flow' - [Bturrock 1971 



Arons & Scharlemann 1979). A concise but to-the-point 



account of essenti al propertie s of the SCLF models has 
been presented by Arons 200C. 
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Three boundary conditions essential to the electro- 
dynamics above the polar cap are (Muslimov & Tsygan 



1992) 



1) £ ■ B = for the magnetosphere within the closed field 
lines, 

2) $ = at the surface and at the interface between the 
closed magnetosphere and the open field lines, 

3) £|| = at the surface level, 

where £m is the electric field component parallel to the 
magnetic field. 

Last but not least, it is assumed that the outflow is sta- 
tionary and the magnetosphere remains axisymmetric. 

The electric field £ required to bring a charged particle 
into corotation satisfies the following equation 



£ + -((« - lj lt ) x r) x B = 

c 



(12) 



where the inertial-frame dragging effect is included 
([Beskin 199(1 Muslimov & Tsygan |1992|) with u; LT = 



K g {R s /r) 3 Cl, where n g = I/(M s Rf) ■ R g /R s , and R g 
2GM s /c 2 . 

Charge density necessary to support this local £ is 



_L_ _ O B 

Qcorot — ~: V ■ C- ~ 

47T Zttc 



The charge density Qcorot due to SCLF at r 
the Goldreich-Julian charge density: 



(13) 



R s is called 



£>GJ 



n b 

2vrc 



(14) 



As charged particles flow out along the open field lines a 
deviation of the local charge density pi OC ai from the local 
corotation density g CO rot develops. By using now two re- 
lations satisfied in the dipolar structure, B(r) cx r~ 3 and 
g(r) cx r -3 , one obtains a simple formula for local devia- 
tion from the corotation charge density: 



f?local Qcorot 



n b 

2ttc 



1 



(15) 



Accordingly, the accelerating potential drop in SCLF 
reads 



A$|| « A$ pc K g 



(16) 



O.I5/45. This is a 



where A$ pc is given by (g) and K t 
remarkable result obtained by Muslimov & Tsygan 1992 : 
the potential drop due to inertial-frame dragging is signif- 
icantly larger than the drop induced by curvature of mag- 
netic field lines ('flaring', Arons fc Scharlcmann 1979] ), or 
- useless to say - by particle inertia (Michel 1974)). 



5. Electric field structure in SCLF gaps 

The model with frame dragging effects presented in a sim- 
plified form above, does not take into account possible 



feed-back effect due to e -pairs formed via photon ab- 
sorption within open magnetic field lines. Copious pair 
creation occurs in a relatively thin layer - a pair forma- 
tion front (PFF). The creation of pairs leads to the screen- 
ing of the accelerating field £11 within the layer of PFF. A 
detailed picture of this effect would require to follow the 
dynamics of electrons and positrons in a self-consistent 
way. Instead, it is reasonable to assume, that the field is 
shorted out at the height were the first e ± -pair is created 
(hereafter denoted as h c ): £|| — for h > h c . 

The problem of electric field structure in the context 
of SCLF with boundary condition £11 = set at ho = 
(stellar surface ) and at h > h c (PFF) was formulated 



and solved by Harding fc Muslimov 1998a The solution 
is rather lengthy and includes special functions. It is how- 
ever possible to obtain simple but quite accurate analyt- 
ical approximations. As long as the length h c of the ac- 
celerator is of the order of polar cap radius R pc , the ac- 
celerating electric field may be approximated according to 



Dyks fc Rudak 2000a as 



£ 



-1.46 



B 



12 



p3/2 



hi 1~y cos X Gauss, (17) 



where B12 = i? s /10 12 G, P is the spin period in seconds, 
X is the angle between the spin axis and the magnetic 
moment of the rotating star, h is expressed in cm, and 
M = 1.4 M©, R s = 10 6 cm. The magnetic colatitude 
£ = 6/6{rj) is scaled with the half-opening angle of the 
polar magnetic flux tube 9(rj), where rj = 1 + h/R s . The 
magnetic colatitude function fx (£) is a monotonically de- 
creasing function, with /i(0) ~ 1 and /i(l) = 0. 

Vertical structure of the electric field depends (via the 
location of PFF) on radiative processes which induce the 
pair creation: curvature radiation (CR) and inverse Comp- 
ton scattering (ICS) on soft X-ray photons from the stellar 
surface (brief characteristics of these processes is presented 
in the next section). 

An interesting way to avoid immediate screening of the 
electric field by created pairs, relevant in strong magnetic 
fields ;> 0.1B cr it, was proposed by Usov fc Melrose 1995 



who noticed that pairs are created then near the energy 
threshold, most favourably in bound states - as atoms of 
positronium - which then move to high altitudes before 
being ionized. 

6. Radiative processes in pulsar magnetospheres 

Cooling of ultrarelativistic electrons via curvature radia- 
tion (CR) and magnetic inverse Compton scattering (ICS) 
are the most natural ways of producing hard gamma-rays 
capable of inducing cascades of e ± -pairs and secondary 
HE photons. These two processes dominate within two 
distinct ranges of Lorentz factors 7 of primary electrons. 

When 7 10 6 , magnetic inverse Compton scat- 
tering plays a dominant role in braking the electrons 
and it is the main source of hard gamma-ray photons 
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( jBturricr ct al. 1995 ). Energy losses due to resonant ICS 
limit the Lorcntz factors of the particles to a level which 
depends on electric field strength £» , temperature T and 
size of hot polar cap, and magnetic field strength B s 
(Xia ct al. 1985, |Sturner 1995 ). The Lorentz factors can 
then be limited even to ~ 10 3 . This stopping effect be- 
comes more efficient for stronger magnetic fields, and it 
was suggested as an explanation for the observed cutoff 
at ~ lOMeV in the spectrum of B1509-58 ( [Sturner 1995| ). 
ICS was also found to be able (under some circumstances) 
to smear out the monotonic energy distribution of elec- 
trons sig nificantly, with interest ing consequences for the 
cascades ( [Dyks fc Rudak 2000b|) . 

However, in their modern versions the accelerators of 
particles are strong enough to outpower the ICS cooling. 
In consequence, very high Lorentz factors - 7 <^ 10 6 - 
are achieved by electrons, limited by CR. The first de- 
tailed scenario of radiative processes in CR-induced cas- 



cades was presented by (Daugherty & Harding 1982) and 



despite many modifications and additions its basic fea- 
tures remain valid. The model assumes that primary elec- 
trons accelerated to ultrarelativistic energies emit curva- 
ture photons which in turn are absorbed by the magnetic 
field and e ± -pairs are created. These pairs cool off in- 
stantly via synchrotron radiation process. Whenever the 
SR photons are energetic enough they may lead to fur- 
ther creation of pairs, etc. ICS can still be incorporated to 
the models with CR-induced cascades as the process rel- 
evant for e ± -pairs, since typical Lorentz factors of theirs 
do not exceed ~ 10 3 . According to the analytical model 



of Zhang fc Harding 2000a the empirical relations for X- 
ray and gamma-ray luminosities of pulsars (presented in 
Sect. 3) can be reproduced satisfactorily when the ICS in- 
volving e -pairs is included. 

Processes relevant for generation and transfer of high- 
energy photons in pulsar magnetospheres are therefore: 

- curvature radiation, 

- magnetic inverse Compton scattering, 

- magnetic pair creation (7 — ► e ) 

- synchrotron radiation, 

- photon splitting (7^7+ 7,) 

- photon-photon pair creation (7 + 7 — » e*). 

Their basic properties are given below, except for 7 + 
7 — * e^: in practice it is treated exactly as in free 
s pace. This is justified in models of 'thick outer gaps' 
( |Chcng fc Zhang 1999 ) but within the framework of po- 
lar cap models this is not the case, in general. However, 
no handy formula is available for the cross section of this 
process in the limit of high B and standard non-magnetic 
formulae are in use (e.g. Zhang & Qiao 199S). 

To illustrate how these processes contribute to the 
high-energy radiation of a pulsar the numerically cal- 
culated components (due to the first four processes 
in the list) are pr esent ed in Figs. | 4| a nd |5] (after 



field in the Vela pulsar does not exceed 10 13 G and pho- 
ton splitting was neglected as non-competitive to the 
magnetic pair creation. The electric field structure of 
the accelerator used in these calculations is taken after 



Harding fc Muslimov 1998a 



6.1. Curvature radiation 

Relativistic electron of energy 7771c 2 (we take 7> 1) slid- 
ing along the magnetic field line of curvature Q ct will emit 
photons with a continuum energy spectrum peaked at the 
characteristic energy 

f 3 

:,,. = -<■/» -U (18) 



3 r 

— — cn — . 

2 Qci 



The radius g cr for a purely dipolar line attached to the 
outer rim of the polar cap can be approximated not far 
away from the NS surface as g CT w ^R s ■ Ric « 10 8 \/Pcm. 
The cooling rate of that electron is 



- _ _ 2 e 2 7 4 

7cr — „ n 

3 mc 



(19) 



For a monoenergetic injection function of electrons Q(j) oc 
5(7 — 70) and their cooling due solely to CR the electrons 
will assume a single power-law distribution in energy space 
7V 7 (el.) oc 7 -4 for 7 < 70 as long as they stay within the re- 
gion of the cooling. Their escape introduces a natural low- 
energy cutoff 7cutoff in Nj(el.). Therefore, the unabsorbed 
CR energy spectrum f s (e) due to the injected electrons 
has a broken power-law shape, with a high-energy limit set 
by 70 and the break at some energy £break- For e > Ebreak 
the energy spectrum is / £ (e) oc e~ 2 / 3 , and / e (e) oc e +1 / 3 
for e < Ebrcak- Since nonthermal spectra cover ususally 
many decades in energy it is more convenient to use sf e (s) 
for easy comparison of power in different parts of energy 
space (see Fig. 0). Accordingly, ef E (e) oc e +1 / 3 above the 
break, and oc e^/ 3 below the break. 

The cutoff limit 7 cu toff can be found by comparing the 
characteristic cooling time scale t CI = 7/|7cr| with the 
estimated time of escape t esc , which we take as t csc ~ 
Qcr/c. Therefore 



£broak« - A n— w 150 MeV, 
4 r 



(20) 



where ro is the classical electron radius (Rudak & 
Dyks 1999). Note, that photon energy Ebreak at which the 



Dyks fc Rudak 2002a| and |Dyks et al. 200l[) along with 
overlaid data points for the Vela pulsar. The dipolar 



cooling break occurs does not depend on any pulsar pa- 
rameters. 

The spectrum of CR calculated numerically to model 
the Vela pulsar (Dyks et al. 2001) is shown in Fig.^ as dot- 
dashed line. High-energy cutoff due to one-photon mag- 
netic absorption occurs around 10 GeV. Note the impor- 
tance of gamma-ray detectors capable to operate above 
10 GeV for (in)validating the model. The low-energy CR 
spectral break £break is prominent at ~ 40MeV. Below 
^break the power of CR decreases and eventually becomes 
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unimportant at 
nent takes over. 



lMeV where the synchrotron compo- 



6.2. Magnetic pair creation 

Pair creation via magnetic photon absorption (7 + B — > 
e ± + B) is kinematically correct because the magnetic 
field can absorb momentum. To ensure high chances for 
the process to occur it is not enough for a photon propa- 
gating at a pitch angle if) to local B to satisfy the energy 
threshold condition, sinip ■ e > 2mc 2 , but high optical 
thickness t~ 1 b of the magnetosphere is required. In fact the 
condition t 1 b = 1 is used as a criterium for the so called 
death-line for radiopulsars in the P — P diagram. Maximal 
values of sin ip for curvature photons in the dipolar field do 



not exceed - 0.1sin6> pc « 0.0014 P -1 / 2 flSturrock 197l| ) 



and, therefore, both conditions arc difficult to be met for 
long-period pulsars. 

The absorp tion coefficient for the process as described 
in Erbcr 1966| and used to calculate r 7 s reads 



m ' 2 A c B, {X} 



'crit 



(21) 



where a is a fine structure constant, A c is a Compton 
wavelength, B crit = m 2 c 3 /en ~ 4.4 x 10 13 G, B± is the 
component of the magnetic field perpendicular to the pho- 
ton momentum, and \ = hj}. is the Erber pa- 



rameter x- The function T(\) is then approximated as 
T( X ) « 0.46exp(-4//3x), which is valid for X ^ 0.2; for 
0.2 this approximation starts to overestimate 77(e) 



X 



> 



significantly. The function / is a near-threshold correc- 
tion introduced by Daugherty k Harding 1983, important 



in the case of high local B. Electron-positron pairs cre- 
ated through the magnetic absorption emit synchrotron 
photons (SR) and also may get involved in ICS with soft 
photons. 

Taking advatage of the dipolar character of the mag- 
netic field, it straightforward to find analytical approxi- 
mation of the condition r 7 s = 1. A photon created at a 
radial distance r from the neutron star center, with its 
momentum parallel to the local magnetic field line will 
undergo magnetic absorption if its energy exceeds the fol- 
lowing local cutoff energy (after Baring k Harding 2000, 
Baring 200 lj) 



.-,,„„„ * (),-| p 1 / 2 ( __ 



1/2 



J 1 0.LB crit ( r . . 
x max { 1, — - I — I > GeV. 



B pc V Rs 



(22) 



This formula is valid for the magnetic field lines with their 
footpoints at the outer rim of the polar cap, of colatitude 
$pc (eq< i|) i hence the dependence on \/~P . For other colat- 
itudes 9 the formula should be multiplied by pc /9. 



6.3. Synchrotron radiation 

Consider a particle of energy jmc 2 gyrating around a local 
field line at a pitch angle ip. Let 711 denotes the Lorentz 
factor of the reference frame comoving with the center of 
the gyration. As long as 711 ^ 1 it relates to the pitch angle 
tp via svaijj ps TiT • The energy available for synchrotron 
emission at the expense of the particle is 7j_ mc 2 , and 7 = 
7-L 7|| • 

The rate of SR cooling reads 



7sr=-^3 2 7i = 

3 m c c 



.H^f^sin^T 2 . 
3 m c c 



(23) 



In comparison to the CR cooling it is enormous (due to 
much smaller curvature radius). 

Critical photon energy (analogous to (|T§|)) reads 



3 * eB 2 • , 
= — n 7 sinifi. 

2 m c c 



(24) 



For a monoenergetic injection function of particles (e - 
pairs in the context of this review) and their cooling due 
to SR the energy spectrum of SR spreads between a high- 
energy limit £ S r(7o) set by 70 of the injected (created) 
particles, and a low-energy turnover e ct determined by 
the condition j± ~ 1: 



e c t 



3 eB 1 
£sr(7 = 711 ) = 77 ft -r 



2 m c c sin-0 



(25) 



The spectrum assumes a single power-law shape f e (s) oc 
£~ x / 2 (and accordingly - ef E (e) oc e +1 / 2 ) above the 
turnover. Below e ct , the spectrum f e changes it slope, 
asymptotically reaching oc e +2 . It is built up by contri- 
butions from low-energy tails emitted by particles with 
7_L 3> 1, and each low-energy tail cuts off at local gyrofre- 
quency, which in the reference frame comoving with the 



center of gyration is ujb = 



eB 



The spectrum of SR 



calculated with Monte-Carlo method to model the Vela 
pulsar is shown in Fig.[I| as a dashed line. The low-energy 
part of the SR spectrum at e ct seems to be essential for 
connecting the RXTE and the OSSE points. 

Synchrotron as well as cyclotron emission become pro- 
hibited at strong local magnetic fields B > 0.1-B cr it 
(|Usov k Melrose 1995[ [Arendt k Eilek 2001[) - the mag- 



netic photon absorption proceeds here close to the thresh- 
old, and the pairs created do not occupy excited Landau 
levels. In consequence, electromagnetic cascades are ex- 
pected to be notably w eaker than for B < 0.1 -B cr jt. Along 
with photon splitting ( Baring k Harding 2001 ), this may 
explain why do radiopulsars avoid the high-B region in 
the P—P diagram. 

6.4. Magnetic Inverse Compton Scattering 

Consider an electron with a Lorentz factor 7 moving along 
a magnetic field line B and a photon of energy e = emc 2 
moving at angle arccos \x to the field line. In the reference 
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frame comoving with the electron (primed symbols) the 
counterpart of the free-space Compton formula, due to 
energy-momentum conservation appropriate for collisions 
with the electron at the ground Landau level both in the 
initial and final state, reads 



e' s =(l-fi' s 2 ) 1 jl + e '(l-///4) + 



1/2 



(26) 



where e' = ey(l — P^i) flHcrold 1997 ), and symbols with no 
subscript and with the subscript 's' refer to the state be- 
fore the scattering and after the scattering, respectively. A 
longitudinal momentum of the electron in the electron rest 
frame changes due to recoil from zero to (e'fx' — e' s /i^)mc. 

The polarization-averaged relativistic magnetic cross 
section in the Thomson regime may be approximated with 
a nonrelativistic formula (Dcrmcr 1990): 



9i 



92 - 9i 



(27) 



where a T is the Thomson cross section, and g\ and 172 are 
given by 



91 (u) = 



92(u) 



(« - 1) 2 + . 



where u = e'/e B , a = 2ae B /3, e B = heB /m 2 c 3 and a 
is a fine-structure constant. The resonance condition for 
the scattering is therefore the cyclotron resonance e' = 
e B . The factor a represents a 'natural' broadening of the 
resonance due to finite lifetime at the excited Landau level. 

In the Klein-Nishina regime (e' > 1) the relativis- 
tic magnetic cross section for the \fi'\ « 1 case be- 
comes better approximated with the well known Klein- 
Nishina relativistic nonmagnetic total cro ss section a KN 
(iDaughcrty fc Harding 1986j |Dcrmcr 1990|) . 

The rate 1Z of scatterings subject by an electron mov- 
ing across the field of soft photons, measured in the lab 
frame is 



/? M ) 



(29) 



where f2 = dfidcj) is the total solid angle subtended by the 
source of soft photons, /i = cos 6, a is a total cross section 
for the process, and dn^Jds/dVl is the local density of the 
soft photons. 

The properties of the field of soft photons are usually 
simplified by taking drip^/ds/dCl as for the blackbody ra- 
diation. This simplification should be taken with care since 
magnetised atmospheres of neutron stars introduce strong 
anisotropy as well as spectral distortions to the outgoing 
radiation (Pavlov et al. 1994). Effectively it means that 



the ICS effects obtained with this simplification are just 
upper limits to the actual effects. 



To estimate electron cooling rate 7 ICS due to the ICS 
the differential form of (E7h is necessary: 



da 



3(T T 

16tt 



(i- M ' 2 )(i- M ; 2 ) + 



+ J (1 + m' 2 )(1 + a4 2 )(5i+52) 



(30) 



(eg. Daughcrty fc Harding 1989) ) , where dQ' s — d(/)' s df/ s is 
an increment of solid angle into which outgoing photons 
with energy e' s in the electron rest frame are directed. The 
mean electron energy loss rate then reads 



7lCS 



de 



( dn ph \ 



(31) 



where e s = £^7(1 + /?aO 1S the scattered photon energy in 
the lab frame (e.g. Dermer 1990). 

The spectrum of magnetic ICS calculated numerically 
to model the Vela pulsar is shown in FigsJ^ and || with 
thin solid line. The blackbody soft photons originating at 
the stellar surface (dotted line) are upscattered by sec- 
ondary e ± -pairs at the expense of their "longitudinal" 
energy 711 mc 2 , assumed to remain unchanged during the 
burst of synchrotron emission. Without the ICS compo- 

(28) 

nent due to the e -pairs the RXTE data for the Vela 
pulsar would be difficult to reproduce within the model. 
It is worth to note, that the magnetic ICS component due 
to primary electrons (not shown in FigJJ) is energetically 
insignificant comparing to the CR component. 

6.5. Photon splitting 

Photon splitting into two photons in the presence of 
magnetic field B i s a third-or der QED process with no 
energy threshold ( Adler 1971 ). The attenuation cocffi- 
cicnt, after averagin g over the polarization states, reads 
(Warding et al. 1997|) 



1 



IOtt 2 A c 



19 
315 



B sin 

B C r\t 



^ — 1 

e cm , 



(32) 



(provided B does not exceed -B C rit substantially) where 
a is a fine structure constant, e is the photon energy in 
units of mc 2 , #kB is the angle between photon momentum 
vector and the local magnetic field. The process, therefore, 
strongly depends on magnetic field strength B. 

Photon splitting has attracted substantial interest in 
recent years due to the discovery of neutron stars with su- 
percritical magnetic fiel ds flMcrcghctti 199E ) . It has been 
analysed in details by Harding et al. 1997 and incorpo- 



rated in a Monte Carlo code tracing the propagation of 
electromagnetic cascades in the magnetospheres of high-_B 
pulsars. The effect was found to explain satisfactorily the 
unusual cut-off observed in the gamma-ray spectrum of 
B1509-58 (see Sect. 3). Generally, it becomes competitive 
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Fig. 4. The model energy spectrum calculated by Dyks fc Rudak 2002a to reproduce the spectral features of the 
Vela pulsar (P = 89 ms, B s = 6 x 10 12 G). The accelerator is located at ho = 3i? s above the surface (see Sect. 6). 
The broad-band spectrum consists of four components due to: curvature radiation of primary electrons (dot-dashed), 
synchrotron radiation of secondary e ± -pairs (dashed), inverse Compton scattering of surface X-ray photons on the 
e ± -pairs (thin solid) and the blackbody surface emission (dotted). The surface temperature T s = 1.26 x 10 6 K was 
assumed for the neutron star. Total spectrum is given by a thick solid line. Phase-averaged data points for Vela from 
different telescopes are indicated: crosses - multicolor photometry with NTT and HST (Mignani & Caraveo 2001); 
filled squares - RXTE (jStrickman et al. 1999Q ; open circles - OSSE ( |5trickman et al. 1996); rectangles - COMP TEL 
(3ch6nfclder et al. 2000); filled circles plus an upper limit just above 10 GcV - EGRET (Thompson et al. 1997). 
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Fig. 5. Same as Fig. || but with two parameters changed drastically (B s — 0.6 x 10 12 G, ho = 0) in order to obtain 
as much power in the optical range as possible by modifying the ICS component due to e ± -pairs (thin solid line). 
Despite unrealistically favourable conditions to down-scatter some of the blackbody photons, the low-energy tail of 
the ICS component falls short by a factor of ~ 3 below the observed optical flux level (crosses). 
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to the magnetic pair creation for dipolar magnetospheres 
with B s *J 0.3-B cr it (Harding et al. 1997). The degradation 



of photon energy in the course of splitting inhibits also any 
development electromagnetic cascades. In consequence, 
high-£? RPP should not emit coherent radio emission. In- 
deed, there exists a high-B region in the P — P diagram 
void of radiopulsars. Even though three recently discov- 
ered (during The Parkes Multibeam Pulsar survey) high- 



B radiopulsars ( Camilo et al. 2000) are located above the 



formal limiting lin e derived in [Baring fe Harding 1998 
(and elaborated in Baring fc Harding 200l| ), the general 
argument that mag netars are expected to be radio-quiet 
RPP remains valid ( [Zhang fc Harding 2000b| ). 



7. Facing the Vela pulsar 

In order to fit our model to the phase averaged spec- 
trum of Vela over almost 8 decades in energy (see Fig.||) 
we were forced to locate the accelerator at the altitude 
of ho ~ 3i? s , i.e. unrealistically high for any version of 
polar-cap models. The same problem was encountered by 
Daughcrty fc Harding 1996 when modelling the CGRO 
data points of Vela. 

Physical justification for lifting the accelerator up 
to ~ 1 R s might rely on a mechanism proposed by 



Harding fc Muslimov 1998a , They considered a sandwich- 
like structure of two pair formation fronts (PFF), con- 
trolled by either CR or ICS, where the lower PFF would 
be formed due to positrons returning from the upper PFF. 
With reasonable surface temperatures (T s ~ 5 x 10 5 K) it 
is ICS which dominates (comparing to CR) in generating 
high-energy quanta at the energetic expense of upward 
moving electrons and downward moving positrons. Such 
a situation is not symmetrical and not stable, therefore. 
Accordingly, the lower PFF tends to elevate until the CR- 
cooling takes over, which happens at ho ~ 1 R s . Unfortu- 
nately, this mechanism - rather promising in the context 
of required ho ~ 3i? s , even though 3i? s is still notably 
more than ~ 1 Rg — has been invalidated by recent numeri- 



cal calculations: Harding & Muslimov 2002 concluded the 



inability of ICS-controlled upper PFFs to develope lower 
PFFs because of too low multiplicities in such cases. 

The requirement of high-altitude accelerator for Vela 
is likely to be extended to other gamma-ray pulsars. A 
piece of hint for such a requirement comes from a strong 
correlation for known gamma-ray pulsars between the in- 
ferred magnetic field strength and the value of high-energy 
cutoff ( Baring fc Harding 200C| ). [Note, that this empirical 
correlation speaks in favor of high- energy cutoffs in pulsar 
spectra originating due to magnetic absorption of gamma- 
rays - a signature of polar- cap models.] This is shown in 
Fig.|| along with the dependence of e C utoff on the location 
of photon emission point r according to eq.(p2|). 

Last but not least - a real challenge for all contempo- 
rary models for Vela comes from the phase-resolved spec- 
tral analysis of pulsed hard X-rays and the phasing of 




B (Gauss) 

Fig. 6. Spectral cutoff energies estimated for eight 
gamma-ray pulsars are plotted against inferred surface 
magnetic field strengths B s . Three shaded bands (note 
the overlaps) indicate the location of e cu toff according to 
eq.(|2^) for three values of the emission altitude (r = 
Rs, 2R S and 3i? s ) and for t he range of spin period s P 
between 0.1 and 5 s. (From Baring fc Harding 2000 and 
Baring 2001; courtesy of Matthew Baring.) 



the lightcurves in optical, X-rays and gamma-rays, pre- 



sented by Strickman et al. 2001 and Harding et al. 2002 



The strong double-peak lightcurve in optical (fig.l in these 
refs.) is particularly enigmatic. The peak-to-peak separa- 
tion is notably smaller than in HE domain. Moreover, the 
leading peak has no counterpart in any other energy range 
(including the radio). The trailing optical peak seems to 
be associated with the trailing RXTE peak and with the 
trailing gamma-ray peak, though the latter lags behind 
the trailing optical by about 0.1 in phase. Extrapolation 
of the gamma-ray spectrum taken for cither of the two 
gamma-ray peaks falls below the phase-averaged optical 
flux level (by a factor ~ 10 in the case of the trailing 
gamma-ray peak). 

The problem of optical emission is also difficult to solve 
in the framework of our phase-averaged approach to model 
the spectrum of Vela. A difficult task we have faced was 
to find a component (without referring to any ingredients 
external to the model) that would account for the optical 
emission, at least in some particular circumstances. FigJ| 
gives an example of what can be achieved at best in the 
optical domain of the model spectrum, but at the price 
paid at the high-energy end of the spectrum. 

8. Millisecond pulsars 

The question of whether millisecond pulsars (i.e. with low 
magnetic fields) do emit gamma-rays will be hopefuly an- 
swered within the next five years with GLAST. On theo- 
retical side, there had been some attempts to address this 
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problem. For instance, Wei et al. 1996 presented broad- 



band spectra of both pulsed and unpulscd emission ex- 
pected in their version of the outer-gap model. In the 
framework of polar-cap scenarios, ^turner fc Dernier 1994 



made their predictions for gamma-ray luminosity £ 7 by 
scaling down with B their analytical formulae for L 7 for 



Rudak & Dyks 1998 


proposed a mod- 


model by 


Daugherty & Harding 1982; 



this modification incorporates a contribution of electron- 
positron pairs to L 7 , which becomes a non-monotonic 



function of B (see Figs. 7 and 8 of Dyks 1998 ). One of 
the consequences of this modification is that expected 
values of L 1 for millisecond pulsars, including J0437- 
4715, are typically about one or two orders of magnitude 
lower than for classical pulsars (see Fig.|| in section |J). 
Quantitatively similar estimate, though for a different rea- 
son, comes from the model of Dermer & Sturncr 1994, If 
these estimates are correct then some millisecond pul- 
sars, including J0437-4715 (see Fig.fj]), should be detected 
with next-generation experiments, providing thus a test- 
ing ground for magnetospheric models relevant for mil- 
lisecond pulsars. GLAST is expected to have sensitivity 
above lOOMeV about 30 times higher then EGRET, and 
it will reach high-energy limit at 300GcV - closing thus for 
the first time the energy gap between the VHE (very high 
energy) domain accessible with ground-based Cherenkov 
techniques and the HE (high energy) domain of satellite 
experiments ( Kamae et al. 2000] ) . Moreover, millisecond 
pulsars are expected to radiate mostly in the range be- 
tween lOGcV and lOOGeV, this class of objects should 
be of interest for t he upcoming IAC telescope MAGIC 
(|Blanch et al. 1998|) . 

For the time being, the only millisecond pulsar 
marginally detected above 100 MeV (at 3.5 a level) 
is PSR J0218+4232, with the spin period P = 2.3 
ms and the inferred m agnetic field B pc ~ 8.6 x 10 8 
G, ( Vcrbunt et al. 1996 ). The inferred luminosity of the 
pulsed emission for 1 steradian opening angle reaches 
L 1 ~ 1.64 x 10 34 erg s _1 ~ 0.07 L sd . Broad-band spec- 
trum of this pulsar differs from the HE spectra known for 
young gamma-ray pulsars: above 100 MeV the photon in- 
dex a p h ~ —2.6 (Kuiper et al. 2000) and the spectrum 
resembles soft spectra of EGRET UID sources; within 
the BeppoSAX range the spectrum is extremely hard, 
with a ph ^ -0.61 ± 0.32 ( |Mineo et al. 2000| ). Off-beam 
geometry (i.e. the line of sight misses the hollow cone 
beam) was recently proposed ( Dyks fc Rudak 2002| ) for 
this pulsar. In such a situation the high energy "cutoff" 
occurs at a relatively low energy ~ 100 MeV (marked by 
two EGRET points, see Fig.||). Unlike in on-beam cases, 
this cutoff is not due to magnetic absorption; it corre- 
sponds to a maximum energy of those electrons which 
emit curvature photons towards the observer. The shape 
of th e apparent cut-off can easily mimic very soft spec- 
tra ( Harding fc Zhang 2001 used similar viewing geom- 




Fig. 7. Cummulative spectral flux of photons expected for 
the millisecond pulsar J0437-4715 at a distance of 140 pc. 
The shaded region shows the range of flux levels due to un- 
certainity in the maximal energy of primary electrons. The 
main part of the spectrum is due to curvature radiation 
of the electrons. The additional feature reaching the VHE 
domain is due to inverse Compton scattering of soft pho- 
tons from the surface with the temperature 4 x 10 5 K. Sen- 
sitivities of EGRET as well as three major HE and VHE 
experiments of the future are also indicated. MAGIC-LZA 
denotes sensi tivity of MAGIC in its Large Zenith Angle 
mode. (After |Bulik et al. 2000| ). 
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etry arguments to explain very soft spectra of EGRET 



-6 -5 -4 -3 -2 -1 1 2 3 4 5 6 
log e [MeV] 

Fig. 8. Two-component theoretical spectrum for the off- 
beam geometry (continuous lines) is compared with the 
spectrum observed for PSR J0218+4232 (diamonds - Bep- 
poSAX, triangles - OSSE, squares - COMPTEL, circles 
- EGRET). The low energy component is due to cur- 
vature radiation from secondary e pairs (with energy 
~ 10 5 MeV acquired upon their acceleration in a resid- 
ual electric field), whereas the high-energy component is 
due to curvature radiation from primary electrons (with 
energy ~ 10 7 MeV). (From |Dyks fc Rudak 2002| .) 



12 



B. Rudak, J. Dyks & T. Bulik: High-energy radiation from pulsars: a challenge to polar-cap models 



UID sources). The BeppoSAX points are explained in 
this model as due to curvature radiation of secondary 
pairs accelerated within the polar gap. Upon the acceler- 
ation in a residual electric field, the e^ 1 pairs should reach 
the energy 711 between 1.5 x 10 5 and 6 x 10 5 , which is 
about 1% of the primary electron energy. A few hundred 
of pairs per primary electron are required to be accelerated 
in order to reproduce the BeppoSAX level. According to 



Harding & Muslimov 2001 and 2002, acceleration of such 



a large number of pairs within the polar gap is difficult to 
achieve. 



9. Viewing geometry effects 

Viewing geometry has been recognized as an important 
factor in shaping the lightcurves as well as broadband en- 
ergy spectra of particular objects (like e.g. the Vela pulsar, 
Harding fc Muslimov 199g| ) as well as in the context 



see 



of EGRET UID galactic sources (e.g. recent review by 
Baring 200rj). 



To taste a variety of apparent pulsar characteristics 
accessible upon changing the viewing angle £ we present 
here two cases: of a classical pulsar in Fig.^ and a millisec- 
ond pulsar in Fig.[l0|. The lightcurves (these are energy- 
integrated above 100 MeV) and the gamma-ray spectra 
were calculated with our numerical code. The spectra are 
two-component: due to curvature and synchrotron emis- 
sion; inverse Compton scattering was neglected. The ef- 
fects of aberration and light travel delays were included. 
These effects are important for millisecond pulsars, and 
affect their lightcurves, especially at large inclination an- 
gles a (the case not shown here). Whenever the line of 
sight misses outer rim of the polar cap, the observed spec- 
tra become soft with an apparent exponential (and not 
super-exponential) decay, mimicking thus an outer-gap 
candidate. The millisecond gamma-ray pulsar J0218+4232 
(Fig.||) is probably a good example of such a case. 

10. Concluding remarks 

High-energy astrophysics of pulsars was challenged by 
unexpected richness of spectral and temporal properties 
found for the brightest gamma-ray pulsars. Numerous 
modifications (both, minor and major) to the existing 
models of magnetospheric activity are being invented to 
accommodate at least some of these properties. It is clear, 
however, that we actually need good quality high-energy 
data for much weaker sources. Only then it will be possi- 
ble to asses on statistical grounds the significance of those 
properties (the argum ent used by m any authors, recently 
reiterated strongly by Baring 2001 ). Chances are that it 
may have happen in a few years time. 



The planned observatory GLAST (Kamae et al. 2000) 




Fig. 9. The characteristics of high energy radiation ex- 
pected from a classical pulsar with P = OAs, B = 10 12 G 
and the magnetic dipole inclination a = 8°. The top left 
panel shows the intensity distribution for radiation above 
100 MeV as a function of rotational phase (f> and viewing 
angle (. Horizontal color-coded lines correspond to the 
position of nine observers located at different angles C- 
For each £ the phase-averaged spectrum drawn with color- 
coded line is shown in top right panel. The corresponding 
pulse profiles are shown in small panels labelled with (. 
The spectrum of total luminosity from the pulsar (thin 
solid line in top right panel) is plotted (at an arbitrary 
level relative to the phas e-averaged spectra) for compari- 
son of the shapes. (From Wozna ct al. 200*2] ). 



will be superior to EGRET in two aspects. First, its sen- 
sitivity at 10 GeV will be more than two orders of magni- 
tude better than that of EGRET. Second, it will reach 



energy of 300 GeV, closing thus for the first time the 
energy gap between satelite and ground-based observa- 
tories ( HE and VHE, res pectively). The MAGIC Tele- 
scope (Blanch et al. 1998) - a 17 m diameter Imaging 
Air Cherenkov Telescope (IACT) - is expected to operate 
with sensitivity about three orders of magnitude higher at 
10 GeV than EGRET. Its advanced technology will make 
possible to cover energy range between 10 GeV and 1 TeV, 
and to reach ~ 50 TeV in the Large Zenith Angle mode. 
Energy ranges of GLAST and MAGIC will overlap over 
more than one decade in energy. 

We expect that the question of whether millisecond 
pulsars emit gamma-rays at all, will be answered relatively 
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advantage of incorporating multipolar components to the 
polar-cap models is twofold: 

1) magnetic field strength may locally be lowered slightly, 
opening thus windows for GeV photons to escape without 
magnetic attenuation; postulating high-altitude accelera- 
tors would be then unnecessary. 

2) directional characteristics of high-energy radiation 
would change notably with respect to the dipolar axis at 
higher altitudes, where radio emission is generated. 
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